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ABSTRACT 

In this work, we investigate the accuracy of various approximate expressions for 
the transit duration of a detached binary against the exact solution, found through 
solving a quartic equation. Additionally, a new concise approximation is derived, which 
offers more accurate results than those currently in the literature. Numerical simula- 
tions are performed to test the accuracy of the various expressions. We find that our 
proposed expression show yields a > 200% improvement in accuracy relative to the 
most previously employed expression. 

We derive a new set of equations for retrieving the lightcurve parameters and 
consider the effect of falsely using circular expressions for eccentric orbits, with par- 
ticularly important consequences for transit surveys. The new expression also allows 
us to propose a new lightcurve fitting parameter set, which minimizes the mutual 
correlations and thus improves computational efficiency. The equation is also readily 
differentiated to provide analytic expressions for the transit duration variation (TDV) 
due to secular variations in the system parameters, for example due to apsidal preces- 
sion induced by perturbing planets. 

Key words: techniques: photometric — planetary systems — eclipses — methods: 
analytical — celestial mechanics 



1 INTRODUCTION 

Transiting exoplanets and eclipsing binaries produce famil- 
iar U- and V-shaped lightcurves with several defining quan- 
tities, such as the mid-eclipse time, eclipse depth and dura- 
tion. Out of these, the transit duration is undoubtedly the 
most difficult observeable to express in terms of the physical 
parameters of the system. Kipping (2008) showed that the 
duration is found by solving a quartic equation, to which ex- 
ists a well-known solution. In general, two roots correspond 
to the primary eclipse and two to the secondary but this 
correspondence determination has an intricate dependency 
on the input parameters for which there currently exists no 
proposed rules. As a consequence, there currently exists no 
single exact expression for the transit duration. 

In many applications, the process of discarding un- 
wanted roots may be performed by a computer, but nat- 
urally this can only be accomplished for case-by-case exam- 
ples. The benefits of a concise, accurate and general expres- 
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sion for the transit duration, as we will refer to it from now 
on, are manifold. The solution provides lower computation 
times, deeper insight into the functional dependence of the 
duration and a decorrelated parameter set for fitting eclipse 
lightcurves (see §6). Such a solution may also be readily dif- 
ferentiated to investigate the effects of secular and periodic 
changes in the system (see §7). With the changes in transit 
duration recently being proposed as a method of detecting 
additional exoplanets in the system (Miralda-Escude 2002; 
Heyl & Gladman 2007) and companion exomoons (Kipping 
2009a,b), there is a strong motivation to ensure an accurate, 
elegant equation is available. 

In this work, we will first propose a new approximate ex- 
pression for the transit duration in §2. In §3, we will derive 
two new approximate expressions for the transit duration 
and discuss others found in the literature; exploring their 
respective physical assumptions and derivations. In §4, we 
present the results of numerical tests of the various formu- 
lae. We find that one of our new proposed expressions offers 
the greatest accuracy out of the candidates. In §5 and §6, 
we utilize the new equation to derive mappings between the 
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Figure 1. Various definitions of the transit duration are marked 
on a model transit lightcurve. 



observeable transit durations and the physical model param- 
eters. These mappings may be used to obtain a decorrelated 
transit lightcurve fitting parameter set, which enhances the 
compuational efficiency of Monte Carlo based methods. Fi- 
nally, in §7, we differentiate the favoured solution to predict 
how the transit duration will change due to apsidal preces- 
sion, nodal precession, in-fall and eccentricity variation. 



2 THE TRANSIT DURATION EQUATION 
2.1 Exact solution 

Before we begin our investigation, let us first clearly define 
what we mean by the transit duration. There exists several 
different definitions of what constitutes as the transit du- 
ration in the exoplanet literature. In figure 1, we present a 
visual comparison of these definitions. There exists at least 7 
different duration definitions: It, £ f, T, W, r, ti2 and £34. tr 
and tp are the definitions used bv lSeager fe Mallen-Ornelad 
l|2003l ) and represent the 1 st to 4 th and 2 nd to 3 rd contact 
durations respectively. These can be understood to be the 
total transit duration and the flat-bottomed transit duration 
(sometimes confusingly referred to as the full duration). T is 
the time for a planet to move from its sky-projected centre 
overlapping with the stellar limb to exitting under the same 
condition. W is a parameter we define in this work as the 
average of tr and £f which we call the transit width. We 
note that T > W since when the planet's centre crosses the 
limb of star, less than one half of the planet's projected sur- 
face blocks the light from the star. However, many sources 
in literature make the approximation T ~ W , which is only 
true for a trapezoid approximated lightcurve (a detailed dis- 
cussion on this is given in §6.3). Finally, t\<z and £34 are the 
ingress and egress durations respectively which are often ap- 
proximated to be equal to one single parameter, r. 

Throughout this paper, we assume that the planet is 
a black sphere emitting no flux crossing the disc of a per- 
fectly spheri cal star. The consequ ences for oblate planets is 
discussed bv lSeager fc"Hul ((2002) and for hot planets with 



significant nightside fluxes in iKipping fc Tinettil (|201Ch . We 
will employ the definition of T for the transit duration. Once 
the equation for T is known, it is trivial to transmute it to 
give any of the other definitions provided in figure 1. 

An exact solution for T, in terms of the true anomaly 
/, is given by integrating d//di between and f a (where 
we use / to describe true anomaly throughout this paper). 
Details of the derivation can be found in K08, but to sum- 
marize we have: 



nf b ,f«) = {^ 7 =L^)-[m)-D(u)} a) 

D(f) = 2 V /T^tan- 1 k/^tan |l - <]~^f 
w; LV l + e 2 J I + ecos/ 



(2) 



D(fb) - D{f a ) = 2\A-e 2 tan- 



Vl — e 2 sin fn 



cos f H + e cos / A 



2e(f - e 2 )sin/ H (ecos/ H + cos/ Af ) 
(1 - e 2 ) sin 2 /m + (e cos fn + cos /m) 



(3) 



Where P is the planetary orbital period, e is the orbital 
eccentricity, D(f) is the 'duration function', /m = (/& + 
fa)/2 and fn = (ft — ,fa)/2. It is clear that two principal 
terms define D and hence we dub the expressions used here 
as the 'two-term' transit duration equation. At this point, 
the outstanding problem is solving for f a and /;,, which are 
the true anomalies at the contact points. 

The K08 solution is derived using Cartesian coordi- 
nates, but we consider here a simpler formulation in terms 
of cos(/). In order to solve, we must first consider the ge- 
ometry of the system. As with K08, an ellipse is rotated 
for argument of pericentre and then orbital inclination. We 
choose to define our coordinate system with the star at the 
origin with the +Z axis pointing at the observer. We also 
choose to align the +X axis towards the ascending node of 
the planetary orbit, which ensures that the longitude of the 
ascending node satisfies = and the longitude of pericen- 
tre (of the planetary orbit) is equal to the argument of the 
pericentre (w = to). In such a coordinate system, we may 
write the sky-projected planet-star separation (in units of 
stellar radii): 



S(f) 



a R 



l + e cos / 



sin 2 isin 2 (o; + /) (4) 



Where aR is the semi- major axis in units of stellar radii 
(a/R t ) and i is the orbital inclination. The two contact 
points occur when S is equal to unity. For a body which 
undergoes both primary and secondary transit, there must 
be at least four solutions, which already is an indication of 
a quartic. If we let c = cos / and rearrange (4) in terms of 
purely terms in c, we obtain the following quartic equation: 
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2.3 Two- Term Expression 



=Qo + Qic + Q 2 c 2 + Q A c + Q 4 c 4 
Qo =(csc 2 i + A 2 (cos 2 lo — esc 2 i)) 2 
Qi =4ecsc 2 i(csc 2 +A 2 (cos 2 uj — cos 2 i)) 
A 2 esc 2 i(e 2 + (2A 2 + e 2 - 2) cos %j) 



Q 



+ 2e 2 csc 4 i(3-A 2 )-2A 4 



Q3 =4e esc 2 i(e 2 esc 2 i — A 2 cos 2u>) 
Q4 = (A 2 cos 2u; — e 2 esc 2 i) 2 + A 4 sin 2 2lo 
A=a R (l~e 2 ) 



(5) 
(6) 
(7) 

(8) 
(9) 
(10) 

(11) 



Equation (5) is a quartic equation not satisfying any 
of the special case quartics which are most easily solved 
(e.g. a biquadratic). Since we have four roots for c = cos/, 
there are eight roots in total for /, of which only four are 
physical. Therefore it is preferable to always work with cos / 
since D(f) may be easily expressed in terms of c as well. 
Although the solutions of a quartic equation are well known, 
two of the roots correspond to the primary transit and two 
to the secondary. The correspondence of which roots relate 
to which contact points varies with an intricate dependency 
on the input parameters. Unfortunately, no known rules or 
relations currently exist for this correspondence and we were 
unable to find a system either. Consequently, there currently 
exists no single exact equation for the transit duration. 

We note that the problem of finding the duration of 
an eclipse is not a new one. An equivalent problem is con- 
sidered in iKopall (| 19591 ) for the time taken for a body to 
move between primary and secondary transit. IKopall (|l959l ') 
showed that an closed-form expression is possible by assum- 
ing i — 90°. However, this assumption would be too erro- 
neous for the purposes of finding the duration of a transit 
event. 



2.2 Approximation Solutions 

In order to avoid the quartic equation, we must make an 
approximation. A useful approximation we can make is that 
£>(/) ~ g c — g(f = / c ), i.e. the planet-star separation is ap- 
proximately a constant value given by the planet-star sepa- 
ration at mid-transit. Defining the transit impact parameter 
as b — a,RQ c cos i, it may be shown that the difference be- 
tween fb and f a is given by: 



sinf H 



ft, - fa \ v / l _ 6 2 ~ 



clrQc smi 



(12) 



In addition to (12), we require (f a + fb)/2. A good ap- 
proximation would appear to be that (fb + fa)/2 ~ f c , which 
is the true anomaly at mid-transit. f c is defined as the point 
where S is a minimum. Differentiating S with respect to / 
and solving for / leads to a quartic expression again and 
thus an exact concise solution remains elusive but a good 
approximation is given by: 



/m — fc 



TV 

2 



(13) 



By combining equations (1) & (3) with (12) & (13), we are 
able to obtain a final expression for the duration, which we 
dub T2 ('two-term'). Testing T2 for the exact solutions for 
/a/ and fn provided precisely the correct transit duration 
for all e, as expected. However, we found that using approx- 
imate entries for these terms severely limited the precision 
of the derived equation for large e (the results of numerical 
tests will be shown later in §4). 

The source of the problem comes from equation (3) 
which consists of taking the difference between two terms. 
Both terms are of comparable magnitude for large e and 
thus we are obtaining a small term by taking the difference 
between two large terms. These kinds of expressions are very 
sensitive to slight errors. In our case, the error is from us- 
ing approximate entries for /m and fn- In this next section, 
we will consider possible 'one-term expressions' which avoid 
the problem of taking the difference of two comparable- 
magnitude terms. 



2.4 One- Term Expression 

There are numerous possible methods for finding 'one-term' 
transit duration expressions. Starting from equation (1), we 
could consider using the same assumption which we used 
to derive the approximate true anomalistic duration, A/; 
i.e. the planet-star separation does not change during the 
transit event. This would yield: 



Ti 



Ti = 



2tt 

p A 



-AS 



TV y/T 



o? r q\ cos 2 i 
anQ c sini 



(14) 
(15) 



Where we have used equation (12) for A/. Another 
derivation would be to assume the planet takes a tangen- 
tial orbital velocity and constant orbital separation from the 
planet, sweeping out an arc of length r c Af. It is trivial to 
show that this argument will lead to precisely the same ex- 
pression for Ti. 



3 CURRENT EXPRESSIONS FOR THE 
TRANSIT DURATION 

3.1 ISeager Mallen-Ornelasl <|2003) Equation 

For a circular orbit, the task of finding the transit duration 
is greatly simplified due to the inherent symmetry of the 
problem and an exact, concise solution is po ssible, as first 
presented bv lSeager fc Mallen-Ornelas! (|2003l ) (SMO03). 



Tsmoo3 = — arcsin 

7T 



clr sin 1 



(16) 



The physical origin of this expression can be seen as 
simply multiplying the reciprocal of the planet's tangential 
orbital velocity (which is a constant for circular orbits), by 
the distance covered over the swept-out arc, oA/. 
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Tsmoo3 = v ■ Ad 



\2iraJ 



(aAf) 



(17) 



Where we expand the first term as the orbital period 
divided by the orbital circumference, and the second term 
as the arc length. It may be shown that: 



A/(e = 0) = 2arcsin 



'R 1 



an sin i 



(18) 



It can be seen that our approximate expression for Af, 
presented in equation (12) is equivalent to equation (18) for 
circular orbits. It is also worth noting that both Ti and T2 
can be shown to reduce down to Tsmoo3 for circular orbits. 



3.2 iTinglev &; Sackettl ([2005) Equation 

iTinelev fc Sacked (|2005h (TS05) presented expressions for 
the duration of an eccentric transiting planet, which has 
been used by numerous authors since (e.g. Ford et al. 2008; 
Jordan & Bakos 2008). It is also forms the basis of a 
lightcurve parameter fitting set proposed by Bakos et al. 
(2007). There are two critical assumptions made in the 
derivation of the TS05 formula. The first of these is that: 

• The planet-star separation, r, is constant during the 
planetary transit event and equals r c 

This is the same assumption made in the derivation 
of the Ti equation. Under this assumption, TS05 quote the 
following expression for T (changing to consistent notation). 



Ttsos = 



r c A<p 



(19) 



Where TS05 define r c as the planet-star separation at 
the moment of mid-transit, v c as the planet's orbital velocity 
at the moment of mid-transit and A<f> as 'the eccentric angle 
between the first and last contacts of transit'. In the stan- 
dard notation, there is no such parameter defined strictly 
as the 'eccentric angle' and thus we initially assumed that 
TS05 were referring to the eccentric anomaly. However, sub- 
stiuting the revevant terms for r c and v c gives: 



2tt y/T 



-.A4> 



(20) 



By comparing (20) to equation (14), it is clear A(j> — Af 
(also note equation (14) was derived under precisely the 
same assumptions as that assumed by TS05 at this stage of 
the derivation). We therefore conclude that the term TS05 
refer to as 'eccentric angle' infact refers to true anomaly. 
This is an important point to make because the derivation 
of the TS05 equation would otherwise be very difficult to 
understand by those working outside of the field. Continue- 
ing the derivation from this point, the second assumption 
made by TS05 is: 

• The planet-star separation is much greater than the 
stellar radius, r>ij, 



Critically, this assumption was not made in the deriva- 
tion of Ti or T2. Using this assumption, TS05 propose that 
(replacing A(f> — > A f to remain consistent with the nota- 
tions used in this work and replacing (1 +p) — > 1 to refer to 
T rather than the duration from contact points 1 to 4): 



A/tsos = arcsin 2 



- a R g c cos^ 1 

CIRQc 



A/t 



a R Q c cos 2 i 



cirQc 



(21) 
(22) 



Where TS05 use equation (22) rather than (21) in the 
final version of T. Therefore, TS05 effectively make a small- 
angle approximation for Af, which is a knock-on effect of 
assuming r 2> R* ■ We argue here that losing the arcsin 
function does not offer any great simplification of the transit 
duration equation but does lead to an unneccessary source 
of error in the resultant expression, in particular for close-in 
orbits, which is common for transits. We also note that even 
equation (21) exhibits differences to equation (12). 

Firstly, inside the arcsin function, the factor of csci is 
missing which is present in both the derivation we presented 
in equation (12) and the derivation of SMO03 for circular 
orbits, equation (18). The absence of this term can be under- 
stood as a result of the r R* assumption. As r — > 00, in 
order to maintain a transit event, we must have i — > (tt/2). 

Secondly, the expression we presented for Af earlier in 
(12) has the factor of 2 present outside of the arcsin function, 
whereas TS05 have this factor inside the function. Further- 
more, the SMO03 derivation also predicts that the factor of 
2 should be outside of the arcsin function and this expres- 
sion is known to be an exact solution for circular orbits. In 
a small angle approximation, arcsin 2x ~ 2 arcsin a;, but we 
point out that moving the factor of 2 to within the arcsin 
function seems to serve no purpose except to invite further 
error into the expression. As a result of these differences, the 
TS05 expression for T does not reduce down to the original 
SMO03 equation and is given by: 



ITS05 



a R g% cos 2 



7T yfT- 



UR 



(23) 



3.3 IWinnl (|2010h Equation 

IWinnl ||201Ch (W10) proposed an expression for T based on 
modification to the SMO03 equation. The first change was to 
modify the impact parameter from a R cos 2 — > q c cir cos i, i.e. 
to allow for the altered planet-star separation for eccentric 
orbits. Secondly, the altered planetary velocity should also 
be incorporated. W10 propose that a reasonable approxima- 
tion for the transit duration is obtained by multiplying the 
SMO03 expressions by the following ratio: 



^(/c)[e = 0] 



9c 



(24) 



Where X is given in IWinnl (|2010h and f c is the true 
anomaly at the centre of the transit. This yields a new tran- 
sit duration equation of: 
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Twin — 



7T VI - e 2 



- afjg 2 cos 2 in 1 / 2 



(25) 



Firstly, we note an obvious improvement of the W10 
expression is that we recover the original SMO03 equation 
for e = 0. Secondly, comparison to the equation for Ti re- 
veals that the two expressions are very similar except for 
the position of an extra g c term. Indeed, the Ti and Twio 
expressions are equivalent in the small-angle approximation. 



4 NUMERICAL INVESTIGATIONS 
4.1 Example Systems 

Insights into the robustness and accuracy of the various ex- 
pressions may be obtained through numerical tests of the 
various approximate expressions. We here compare the ac- 
curacy of four expressions: Trsos, Twio, Ti and Ti. These ex- 
pressions depend only on five parameters P, a/R,, b, esinu 
and ecosw. One of the clearest ways of visually comparing 
the equations is to consider a typical transiting exoplanet 
example with system parameters for a/ R* and b and vary 
the eccentricity parameters. P may be selected by simply 
assuming a star of Solar density. 

We create a 1000 by 1000 grid of e sin cj and e cos uu val- 
ues from -1 to 1 in equal steps. Grid positions for hyperbolic 
orbits (e > 1) are excluded. We then calculate the transit 
duration through the exact solution of the quartic equation, 
Tko8, plus all four approximate formulae. We then calcu- 
late the fractional deviation of each equation from the true 
solution using: 



V, 



<- candidate 



Tko8 



candidate 



(26) 



We then plot the loci of points for which the deviation 
is less than 1% (i.e. 2? candidatc < 0.01). In figure 2, we show 
four such plots for different choices of a/R* and b. The plot 
reveals several interesting features: 

• Ti consistently yields the largest loci. 

• Tb is sometimes accurate and sometimes not, support- 
ing the hypothesis that the approximation is not stable. 

• Twio also yields consistently large loci. 

• Ttso5 consistently yields the smallest loci. 



define the 'improvement' of the Ti expression relative to the 
Ttso5 equation as: 



Tri = [(£>TS05/2>ri) - 1] * 100 



(27) 



Where X is measured in %. We can see that if TS05 gives 
a lower deviation (i.e. more accurate solution), we will obtain 
X < 0% whereas if the candidate expression gives a closer 
solution we obtain X > 0% and is essentially the percentage 
improvement in accuracy. For the range < e < 1, we find 
the mean value of this parameter is Tti = 210% and for the 
range < 0.9 < e < 1.0 we find X T i = 458%. We note that 
one caveat of these tests is that they are sensitive to the a 
priori inputs. 

For the case o f Kepl er photometry, following the method 
of lKipping et alj (|2009d ). we estimate that the typical mea- 
surement uncertainty on T will be ~0.1% in most cases. We 
find that Ti is accurate to 0.1% or better over a range of 
|esinw| < 0.5 and |ecosu>| < 0.85 on average. 



5 ANALYTIC INVESTIGATIONS 

5.1 The Consequences of Using Circular 
Expressions for Eccentric Orbits 

SMO03 showed that the 1 st to 4 th contact duration, t T , and 
the 2 nd the 3 rd contact duration, ij?, may be used to derive 
an, i, b and the stellar density, p*. We here consider how 
biased these retrieved parameters would be if we used the 
circular equations for an eccentric orbit. From here, we will 
employ the Ti expression for the transit duration, as this 
equation has been shown to provide the greatest accuracy 
in the previous section. According, to SMO03, the circular 
transit durations are given by: 



t T (SMO03) 
i F (SMO03) 



P ( J(l+p) 2 - a 2 cos 2 ; 

— arcsm — : — : 

n \ an sin i 



P I \/7T— p) 2 — a 2 „ cos 2 i 

— arcsm ; — : — 

n \ an sin i 



(28) 
(29) 



Modification of the Ti solution gives: 



4.2 Additional Tests 

We briefly discuss additional tests we performed for two sets 
of 10 7 different hypothetical transiting exoplanet systems; 
one for eccentricities 0.0 < e < 1.0 and one for 0.9 < e < 1.0. 
In all cases, we randomly generated the system parameters 
weighted by the transit probability and calculated the devi- 
ation of the various formulae. 

We found that the Ti expressions was consistently the 
most accurate, with the W10 of similar accuracy but higher 
assymmetry. We therefore find that the results yield a overall 
preference for the Ti approximation. We note that authors 
using the Twio formulation can also expect an extremely 
good approximation but for the remainder of this paper we 
will only consider using Ti for the later derivations. We may 



p qI . 

it 1 — , arcsin 

7T VI -e 2 

P Qc 
if 1 — , arcsin 



y/Q- +p) 2 - a \ei cos 2 ; 

anQc sini 

y/(l -p) 2 - a 2 R g 2 cos 2 ? 
a R Q c sini 



(30) 
(31) 



Using (28) and (29) , SMO03 show that the impact pa- 
rameter may be retrieved by using: 



[bdcrivod(SMO03)] 2 = 



fl _ n \ 2 - an (tF-n/n , 
I 1 P> sin^(t T 7r/P) V 1 + 

i sin2(t F7r /P) 



(32) 



sin 2 (t T -/r/P) 

Using the same equations for an eccentric orbit gives: 
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m a/R = 5 
b = 0,5 




b = 0,75 






i 




1 







1.0 



-0.5 



a/R = 1Q. 
b = 0,7 



0.0 
esintj 



0.5 



1.0 




Figure 2. Loci of points for which the accuracy is better than 99% for all four candidate expressions, as a function of eccentricity. The 
Ti expression offers both consistency and excellent accuracy. Other system parameters fixed to typical transit values. Blue is for Txso5> 
black is for Twio, re d * s f or ?l an< ^ green for T^. The gray ellipse represents the allowed physical limits. 



[bdc ri vod(SMO03)] 2 = 1 + p 2 + 2p- 



• arcsin( 



V(i- P ) 2 -i> 2 . 

clr g c sin i 



+ sm 



. arcsin( 



V(i+ P ) 2 -ft 2 

a RQc sin i 



angle approximations, this yields ^derived — b 2 . However, for 
larger e sin oj and e cos oj values, the overall effect is to over- 
estimate b for eccentric orbits. 

In addition to the impact parameter, SMO03 proposed 
the parameter cir = a/R, may be derived using: 



V(i-p) 2 -i> 2 - 

»i cir q c sin i > 



= arcsinl — '—. — : — 

,2 a RQc sm 1 ' 

(33) 



derived 

(SMO03)] 2 



(1+P) 2 



Where it is understood that for terms on the right-hand 
side with b in them, we are referring to the true impact pa- 
rameter, b = aRQ c cos i. We plot this function in the case of 
an — 10, b 2 = 0.5 and p = 0.1 in figure 3. Making small- 



sin 2 (t T n/P) 



derived 

(34) 



If we use the assumption bderived 
yields: 



b, then this equation 
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Figure 3. // one uses the circular expressions, the retrieved im- 
pact parameter (squared) is heavily biased by eccentricity. In this 
example, the true value of b 2 is 0. 5 but the introduction of eccen- 
tricity causes b 2 to be overestimated. 




Figure 5. If one uses the circular expressions, the retrieved value 
of p* is heavily biased by eccentricity. In this example, the true 
value of p* is lpQ but the introduction of eccentricity causes p* 
to be underestimated. 



(a/fl i d „ 




Figure 4. // one uses the circular expressions, the retrieved value 
ofa/R* is heavily biased by eccentricity. In this example, the true 
value of a/R t is 10 but the introduction of eccentricity causes 
a/R* to be underestimated. 



[^.derived (SMO03)] 2 = b 2 + 
((l+p 2 )-fe 2 )csc 2 ' 



^/(l+p) 2 -b 2 

angc sini 



(35) 



With small-angle approximations, we have: 



an, derived — a R i / gi cos A i + 



(1 



(36) 



The term inside the square root goes to unity for cir- 
cular orbits, as expected. The deviation in an can be seen 
to become quite significant for eccentric orbits, as seen in 
figure 4 where the exact expression for (34) is plotted. This 
will have significant consequences for our next parameter, 
the stellar density. Stellar density is related to an by ma- 
nipulation of Kepler's Laws: 



p*+p pp 



P* 



3n 3 

Gi^ aR 

37T 3 

GP^ aR 



(37) 



Where the approximation is made using the assumption 
p <C 1. We can therefore see that: 



P*, derived — P* 



2 2 • , (1- e )sin i 
g c cos % H 



P*, derived — — 



(f + esinu) 



3/2 



(1_ e 2)3/2 



(38) 
(39) 



Where in the second line we have assumed that i ~ tt/2. 
A series expansion of ^ into first order of e yields ^ ~ 1 + 
3e sin uj+0[e ] . So observers neglecting an eccentricity of e ~ 
0.1 may alter the stellar density by 30%. As an example, if we 
decreased the density of a solar type G2V star by 30%, the 
biased average stellar density would be more consistent with 
a star of spectral type KOV. Indeed, asteroseismologically 
determined stellar densities of transiting systems could be 
used to infer 

This density bias, which is plotted in figure 5, could be 
extremely crucial in the search for transiting planets. Many 
discovery papers of new transiting planets have only sparse 
radial velocity data and usually no secondary eclipse mea- 
surement. As a result, the uncertainity on the eccentricity 
is very large. 

Critically, planets are often accepted or rejected as be- 
ing genuine or not on the basis of this lightcurve derived 
stellar density. If the lightcurve derived stellar density is 
very different from the combination of stellar evolution and 
spectroscopic determination, these candidates are generally 
regarded as unphysical. This method of discriminating be- 
tween genuine planets an d blends, which may mim i c such 
objects, was proposed by ISeager fe Mallcn ^Qrnelas! (|2003| ) 
(see §6.3 of SMO03) but crucially only for circular orbits. 

Since the typical upper limit on e is around 0.1 in dis- 
covery papers, then the lightcurve derived stellar density 
also has a maximum possible error of ~ 30%. In practice, 
the uncertainty on e will result in a larger uncertainty in p» . 
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Typical procedure is to fix e = if the radial velocity data 
is quite poor, despite the fact the upper limit on e ~ 0.1. As 
a result, the posterior distribution of p» would be artificially 
narrow and erroneous if e 7^ 0. We propose that global fits 
should allow e to vary when analyzing radial velocity and 
transit photometry, as well as a fixed e = fit for compari- 
son. This would allow the full range of possible eccentricities 
to be explored, which would result in a broader and more 
accurate distribution for b, an and critically p*. 



5.2 General Solution for the Physical Parameters 

In the previous subsection we saw how using the circular 
expressions to derive an and p» can lead to severe errors for 
even mildly eccentric systems. We here present expressions 
which will recover excellent approximations values for b, an 
and p„. The new equations are given by: 



(1-P) 2 



sin [(ip7r\ 



-e 2 )/(Pgg)] 



sin 2 [(iy 7T"\ 



2 )/(Pe 2 c )] 



sin 2 [(t F ir-\ 



e 2 )/(Pel)] 



sin 2 [(tj^iry 



2 

a.R = 



(1+P 2 ) 



o • 2 

QZ sin 



37T 3 

T a R 



[(i T7 rvT^?)/(Pf??)] Ql 



GP 2 



P 3 pp 



(40) 



(41) 



(42) 



These expressions can be shown to reduce down to the 
original SMO03 derivations if e — >■ (equations (7) & (8) 
of SMO03). The new stellar density parameter may be used 
with floating e and ui values to correctly estimate the prob- 
ability distribution of this critical parameter. 



W(e = 0) = - arcsin ( ^ 1 



(45) 



ce = an sm 1 



[(l+p) 2 -b 2 ] 1/2 [l- 



(1-p 2 ) + b 2 ] 1/2 - [(1 - p) 2 - b 2 ] 1/2 [l - (1 +p 2 ) + b 2 ] 1 ' 2 



(46) 



It may be shown that a 7^ y/1 — b 2 / (an sini) and thus 
W 7^ T. In the same manner as we derived W , r may also 
be written as a combination of the relevant arcsin functions. 
However, we can already see that such an expression will 
also be extremely laborius. This means that writing down 
the expressions for b(T,r) and an(T,r) is much more chal- 
lenging than that for b(tr,tF) and afl(tT,tf) and we were 
unable to find an exact inversion relation. Indeed, finding 
exact expressions for these inverse relations is unncessary as 
we may retrieve b and an using equations (40) and (41). 

For a circular orbit and in the limit an 2> 1, the relative 
difference betwen T and W may be written as: 



Ti - Wi 2^T _ &2 - ^{l+p) 2 -b 2 - y/(l-p) 2 -b 2 



2^1^62 



(47) 

Notice how for p — > this expression yields zero, which 
is expected since the planet now takes infintessimal size. The 
denominator also reveals that the difference diverges rapidly 
for near-grazing transits i.e. b — > 1. 



5.3 The Transit 'Width' Duration 

In the previous subsection, we have derived the physical pa- 
rameters in terms of tr and £f- We may naively assume 
that this is interchangeable with expressions in terms of T 
and t. Assuming the transit lightcurve is symmetric (exactly 
valid for circular orbits and a very good approximation for 
eccentric orbits), the following relations between these two 
definitions exist: 



tr — tp 



T ^tT+t L=w 



(43) 
(44) 



For the latter, T = W for the trapezoid approximated 
lightcurve only. This is because T is defined as when the 
sky-projection of the planet's centre is touching the stellar 
limb i.e. S = 1. At this point, the fraction of the planetary 
disc occulting the stellar disc is not equal to one half of the 
total in-transit occulted area. In contrast, we here define 
W as the duration between the midway of the ingress to 
the midway of the egress. We can intuitionally see at the 
moment S = 1, less than half of the total area must be 
occulted and therefore T > W = (tr + t F )/2. 

A further validation of this can be seen by writing down 
the equations for tr and tp and combining them using arcsin 
trigometric identities. For the simple case of a circular orbit, 
the resultant expression would give: 



6 APPLICATIONS TO LIGHTCURVE 
FITTING 

6.1 Fitting Parameter Sets: {t c , p 2 , Wi, n} 

The transit lightcurve is essentially described by four phys- 
ical parameters, which form the parameter set {t c , p, an, 
6}. However, efforts to fit transit lightcurves using this pa- 
rameter set is known to be highly inefficient due to large 
i nter-paramete r corre lations, in particular between an and 
b. lCarter et al.l (|2008l ) used a Fisher analysis to show that for 
a symmetric lightcurve, which is approximated as a piece- 
wise linear model (i.e. a trapezoid), a superior parameter 
set to fit for is given by {t c , p 2 , T, r} where r is the 
ingress or equivalently egress duration assuming a symmet- 
ric lightcurve. In our case, these parameters become {t c , p 2 , 
Ti, ti}. The authors reported that using this parameter set 
decreased the correlation lengths in a Markov Chain Monte 
Carlo (MC MC) fit by a factor o f ~ 150. 

In the ICarter et al.l (|2008h analysis, the lightcurve is 
a symmetric trapezoid and therefore W = T. As we have 
already seen, this is not true for a real lightcurve. This raises 
the ambiguity as to whether this fitting parameter should be 
W or T for real lightcurves. 

One advantage of the T parameter is that it is indepen- 
dent of p, whereas W is not. With one degree less of freedom 
than W, T will always exhibit lower correlations and may 
be determined to lower uncertainity. This makes T ideal for 
transit duration variations (TDV) studies, as pointed out 
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bv lKipping et al l |2009d) . However, as we saw in §5.3, there 
presently exists no known expression for converting T and 
t into the physical parameter set which is used to actually 
generate a model lightcurve. 

When fitting a transit lightcurve, our hypothetical al- 
gorithm must make trial guesses for 'the fitting parame- 
ter set' which is then mapped into 'the physical parameter 
set'. These physical parameters are then fed into a tran- 
sit lightcurve model generator, allowing for the goodness- 
of-fit between the trial model and the observations to be 
made. This mapping proceedure is unavoidable since the 
transit lightcurve is essentially generated by feeding the sky- 
projected planet-star separation, 5*, as a f unction of time, 
into a lightcurve generating code like that of lMandel fc Ago] 
l|2002h . Since S (equation 4) is a function of the physical pa- 
rameter set, and not the fitting parameter set, the mapping 
between the two sets is a pre-requisite for any lightcurve 
fitting algorithm. 

Unless an approximation is made that Ti ~ Wi, there 
currently exists no expressions which perform this mapping 
proceedure for the fitting parameter set {t c , p 2 , Ti, ti]q 
Specifically, there currently exists no exact expression for 
&(Ti,ri,p) and clr(Ti,ti,p). Therefore, {t c , p 2 , Ti, n} can- 
not be used as a fitting parameter set unless we assume 
Ti ~ Wi and use {t c , p 2 , Wi, n} 

Fortunately, the consequences of making this assump- 
tion will not be severe, in most cases. This is because the 
trial fitting parameters serve only one function - to produce 
trial physical parameters. These trial physical parameters 
may be slightly offset from the exact mapping but this is 
not particularly crippling since the model lightcurve is still 
generated exactly based upon these trial physical parame- 
ters. The only negative consequence of using this method is 
that an additional correlation has been introduced into the 
fitting algorithm since the offset between T and W will be 
a function of b and p. This correlation will be largest for 
near-grazing transits since equation (47) tends towards oo 
as b — > 1. In general, we wish to avoid such correlations as 
much as possible to allow the algorithm to most efficiently 
explore the parameter space. 



6.2 Fitting Parameter Sets: {t c , G\, Wi, Ai} 

For a trapezoid approximated lightcurve, ICarter et al.l 
(2008) showed that the fitting parameter correlations are 
decreased further by using the parameter set {t c , G, W, 
A}, where A is the area of the trapezoid-approximated 
lightcurve, and G is the gradient of the ingress/egress (note 
we have changed the original notation from S to G to avoid 
confusion with equation 40 

The area of the trapezoid lightcurve is given by S(tr + 
£f)/2 where 8 — p 2 . The gradient of a trapezoid slope is 
given by S/t. Since r = (tr — £f)/2 then both A and G may 
be written as a function of t-r and tp only, thus obviating the 



use of T and the associated issues discussed in the previous 
subsection. 

In order to proceed, a mapping from {G, W, A} — > {p, 
an, b} is required for accomplishing this goal. The exact 
solutions for G, W & A may be found by solving the quartic 
equation discussed in §2. However, the roots of this equation 
yields G(p, an, b), W(p, an, b) & A(p, an, b) whereas we need 
the inverse relations. Since no concise analytic solution for 
the inverse relations currently exists, these inverse relations 
would have to be calculated through a numerical iteration 
but such a process would need to be repeated for every single 
trial leading to vastly greater computation time for a fitting 
algorithm. 

Therefore, a practical compromise is to use approx- 
imate formulae Gi(p,an,b), Wi(p,an,b) & Ai(p,an,b), 
which are easily manipulated to give the inverse relations: 
p(C?x,Wi,Ai), a H (Gi,Wi,Ai) and 6(Gi,Wi,Ai). 

It is critical to understand that using equations for a cir- 
cular orbit or an approximate eccentric expression of poor 
accuracy will cause fitting algorithms to wander into unphys- 
ical solutions and/or increase inter-parameter correlations 
for planets which are eccentric, near-grazing, very close-in, 
etc. It is therefore imperative to use a mapping which is as 
accurate as possible in order to have a robust fitting algo- 
rithm. The mappings to convert the trial {Gi, Wi, Ai} into 
the physical parameters {p, an, b} are given by equations 
(40) and (41) combined with the following replacements: 



/Ai 

p = Vm 

tr,i =Wi + 
t F ,i =Wi- 



Ar 
WxGi 

Aj 
WiGi 



(48) 
(49) 
(50) 



We note that the favoured parameter set derived by 
ICarter et al.l (|200ct ) assumed a symmetric lightcurve which 
is not strictly true for e > 0. However, for an eccentric or- 
bit, the degree of assymetry between the ingress and egress 
is known to be very small (K08, W10) and thus may be 
neglected for the purposes of choosing an ideal fitting pa- 
rameter set. 



6.3 Fitting Parameter Sets: {t c , p 2 , C/R*, b 2 } 

The two parameter sets proposed so far have both included 
W\. Since we know Wi is a function of p but Ti is not, 
any parameter set using W\ will likely exhibit larger corre- 
lations since there is an extra parameter dependency. How- 
ever, a parameter set based upon T\ would have to satisfy 
the condition that it can be inverse mapped into the physical 
parameters. 

A s earch through the li terature finds just such a param- 
eter set. iBakos et all (|2007T) proposed the fitting parameter 
set {t c , p 2 , C/R#, b 2 }, where (,/R* is defined by: 



Although mappings have been proposed, they all make the as- 
sumption T = W 

2 We note that ICarter et al. proposed an additional 

slightly improved parameter-set, but this set required reliable 
prior estimates of b and lacked a physical interpretation 



R, 
R, 



TtS05 

27ra 



1 + e sin uj 



PR* v / i —- ?VT _- 62 



(51) 



C/R; originally defined bv lMurrav fc Dermottl < |l999h . 
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can be seen to be reciprocal of one half of the transit dura- 
tion as defined by TS05. Unlike the {t c , p 2 , T, r} parameter 
set, we do not need to assume T = W to produce an inverse 
mapping. By using Ttsos and b 2 , an exact inverse mapping 
to the physical parameters is possible which offers significant 
advantages. 

Having satisfied the criteria of being both a decorrelated 
parameter and inverseable, C/R* would appear to an excel- 
lent candidate for lightcurve fitting. A further improvement 
to this parameter is possible by using the new Tx approxi- 
mation for the transit duration. Let us define: 



2 



_J_ _ 2tt yT 

~rZ ~ T I 



sj\ — a 2 R g'j. cos 2 i 
a R g c sin i 



The inverse mapping would use the expression: 



2 1 

a R = — 



2 

- — esc 



r 27rVt - e 2 



(52) 



(53) 



Table 1. For each proposed lightcurve fitting parameter set (left 
column), we calculate the inter-parameter correlation matrices 
in the examples of i) a hypothetical near-grazing hot-Jupiter on 
a circular orbit ii) a system similar to the eccentric planet HD 
80606b. We diagnolize the correlation matrices to give M and 
then quantify the departure from a perfectly optimal parameter 
set (right column), where it is understood that corresponds to 
optimal and larger values correspond to greater inter-parameter 
correlations. 



Parameter Set 


Eti - 1| 


Circular orbit example 




{tc p, a R , b} 


2.19333 


{t c ,p 2 , C/R*,b 2 } 


1.71236 


{t c ,p 2 , T/R,, b 2 } 


1.32974 


{t c ,p 2 , Wx, n} 


1.67485 


{t c , Gx, Wx, Ax} 


2.23820 


Eccentric orbit example 




{t c , p, a R , b} 


2.46676 


{t c ,p 2 ,(/R*,b 2 } 


1.56816 


{t c> p 2 , r/R*, b 2 } 


1.56776 


{t c , p 2 , W X , Tx) 


1.57730 


{to, Gx, Wx, Ax} 


2.52948 



6.4 Circular Orbit Example 

Despite the analytic arguments made so far, the clearest 
validation of which fitting parameter set to employ may be 
resolved through numerical simulations. This may be done 
by considering an example system, generating a lightcurve, 
adding noise and then refitting using an MCMC routine 
which outputs the inter-parameter correlations. We note 
that analytic expressions for the covariances may be found 
through Fisher information analysis through the calcula- 
tion of the relevant partial derivatives (|Pall |2008| ). How- 
ever, the equat i ons d escribing the lightcurve, as given by 
IMandel fc Ago! (|2002l ) are quite elaborate and such an anal- 
ysis remains outside of the scope of this paper. Currently, 
there exists no exact Fisher information analysis in the liter- 
ature to draw upon. ICarter et al.l (|200ct ) avoided this prob- 
lem by making a trapezoid approximation of the lightcurve 
and then implementing a Fisher analysis. As discussed ear- 
lier, this requires that we assume T = W, which in itself 
introduces a host of correlations which would be missed by 
the Fisher analysis methodology. Therefore, exact numerical 
testing provides a useful alternative to avoid these issues. 

First, we consider a super-hot Jupiter on a circular or- 
bit with a planet-star separation of a R = 3.5 from a Sun- 
like star (P — 0.76days). We choose to consider a near- 
grazing transit with 6 = 0.9 corresponding to an orbital 
inclination of 7 5. 1°. T he lightcurve is generated using the 
IMandel fc AtoI |2002h algorithm with no limb-darkening 
and 0.25mmag Gaussian noise over a 30 second cadence. 
The lightcurve is then passed onto a MCMC fitting algo- 
rithm where we try several different parameter sets: 

• {to, p, a R , b}: the physical parameter set. 

• {to , p 2 , C/R*, b 2 }: a su ggested lightcurve fitting param- 
eter bv lBakos et al.l (|2007r i. based upon the TS05 duration 
expressions. 

• {t c , p 2 , T/R*, b 2 } : a mo dified form of the fitting pa- 
rameter bv lBakos et al.l |2007), accounting for the improved 
approximate expression for T. 



• {t c , p 2 , Wx, rx}: a suggested set bv lCarter et al.li|200Sl ). 
where Wx and Tx are calculated using the expressions pre- 
sented in this paper. 

• {to, Gx , Wx , Ax}: a second suggested set by 
ICarter et al.l |2008t ). where Gx, Wx and Ax are calculated 
using the expressions presented in this paper. 

In the MCMC runs, we set the jump sizes to be equal 
to ~ 1-er uncertainties from a preliminary short-run. We 
then start the MCMC from 5-<r's away from the solution for 
each parameter, and use 500,000 trials with a 100,000 burn- 
in time. We then compute the cross-correlations between 
the various parameters in trials which are within A^ 2 = 1 
of X 2 1 best (errors rescaled such that \ 2 1 best equals number 
of data points minus the degrees of freedom). We calculate 
the inter-parameter correlations and construct correlation 
matrices for each parameter fitting set. As an example, the 
correlation matrix for the {t c ,p, an, b} parameter set is given 
by: 

CoTi({t c ,p,a R ,b},{to,p, a R ,b}) = 

/ 1 Corr(t. c ,p) Corr(t c ,a R ) Con(t c ,b) \ 

Corr(p,f c ) 1 Corr(p, a R ) Corr(p, b) 

Corr(a H ,i c ) Corr(a fl ,p) 1 Corr(a fl ,fe) 

\ Corr(fe, t c ) Corr(b,p) Corr(p, a R ) 1 J 

We then calculate the semi-principal axes of correla- 
tion ellipsoid by diagnolizing the matrices. For a completely 
optimal parameter set, this diagnolized matrix would be 
the identity matrix. We quantify the departure of each pro- 
posed parameter set from the identity matrix by calculating 
53 -=i — 1| wnere M is the diagnolized correlation ma- 
trix. We display the results in upper half of table 1. 

The correlations of the physical parameter set are pre- 
dictably very large, in particular between a R and b which 
approaches unity. An inspection of the correlations for the 
other proposed parameter sets suggests that the modified 
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Figure 6. Comparison of the correlations between C/-R* against 
b 2 and T/i?* against b 2 . Data comes from fitting a synthetic 
hot-Jupiter lightcurve on a circular, near-grazing orbit with an 
MCMC algorithm. The new T/R* parameter provides two-fold 
lower correlation and preserves the ability to be inversely mapped 
to more physical parameters. The three different types of shading 
represent the 1-a , 2-a and 3-cr confidence regions. 

iBakos et al.1 l|2007l ) formulation offers the lowest correla- 
tions. The effect of modifying C/-R* to T/R* produces a 
clear improvement in the corresponding correlations, as seen 
in figure 6. As a resul t, the numerical tes ts support using the 
modified form of the B akos et alj (|2007T ) parameter set. 

6.5 Eccentric Orbit Example 

As a second example, we consider a highly eccentric orbit. In 
this case, we choose to use the real system HD 80606b. HD 
80606b is highly eccentric planet with e ~ 0.93 first discov- 
ered through radial velocity (Naef et al. 2001) and then later 
found to go through both secondary and primary eclipse 
(Laughlin et al. 2005; Fosse y et al. 20 09). Using the sys- 
tem parameters from I Winn et al.1 ([2009) , we generate a syn- 
thetic lightcurve of cadence 1 minute and 1.0 mmag Gaus- 
sian noise. We ensure the quantity of out-of-transit baseline 
data is approximately equal to the full transit duration. We 
adopt the same methodology as before to calculate the inter- 
parameter correlations of the various sets with the results 
shown in the lower half of table 1. 

In these tests, we find three sets produce approximately 
the same opt imization but the lo west correlations occur for 
the modified IBakos et al.l (|2007T ) set again. The difference 
betw een the corr e lation s in the unmodified and the mod- 
ified IBakos et al.l ()2007l ) parameter set is extremely small, 
but there is a very slight improvement in the modified ver- 
sion. The difference for circular orbits was larger due to the 
more grazing transit and the fact the TS05 expressions do 
not reduce down to the circular form . In co nclusion, we advo- 
cate using the modified lBakos et al.l |2007l ) parameter set i.e. 
{t c ,p 2 , T/Rt, b 2 } to most efficiently fit transit lightcurves. 



7 SECULAR TRANSIT DURATION 
VARIATION 

Transit duration variation (TDV) can occur two possible 
formats: i) periodic change ii) secular change. As an exam- 
ple, periodic change in the transit duration was predicted 
to occur for a transiting planet with a companion moon by 
Kipping (2009a, b). Conceptually analagous to the radial 
velocity of finding planets, the sky-projected tangential ve- 
locity of the planet oscillates around some local mean value 
as a result of the moon's gravitational tug. These changes in 



tangential velocity induce variations in the transit duration, 
allowing for the detection of exomoons as small as 0.2Me 
with space-based photometry (Kipping et al. 2009c). 

Secular changes in tra nsit duration can be ca used by 
numerous possible scenarios. IJordan fc Bakosl (|2008T l showed 
that apsidal precession would induce changes in T and used 
the TS05 equation to predict the size of these changes. As 
we have already demonstrated a better formulation for T is 
possible; we will here present an improved equation for the 
rate of change in T due to apsidal precession, or essentially 
cha nges in uj. 

I Jordan fc Bakosl l)2008h argued that apsidal preces- 
sion can be caused by stellar oblateness, general rela- 
tivistic effects and/or a perturbing planet. Additionally, 
iMurrav fc Dermottl ij 19991 ) showed that in general nodal pre- 
cession should also occur whenever apsidal precession oc- 
curs, leading to changes in the orbital inclination angle, i. 
These changes lead to another form of secular TDV. 

Additionally, we consider here that falling planets, such 
as proposed for WASP-18b (Hellier et al. 2009), would ex- 
perience a changing semi-major axis, a, leading to another 
form of secular TDV. Finally, we will consider the effect of 
varying the orbital eccenticity. All four possible TDVs will 
be derived here using T\ , since this expression demonstrated 
the greatest accuracy in numerical tests. 

7.1 Apsidal Precession 

Apsidal precession is the precession of the argument of peri- 
apse over time and it may be induced from several different 
effects including: 

• General relativistic effects (Einstein 1915; Pal & Kocsis 
2008) 

• Rotational quadrupole bulges on the planet (Sterne 
1939) 

• Tides raised on the planet and the star (Sterne 1939) 

• Stellar quadrupole moment (Murray & Dermott 1999) 

• Kozai mechanism (Kozai 1962) 

• Perturbing planets (Murray & Dermott 1999; Miralda- 
Escude 2002; Heyl & Gladman 2007) 

Out of these examples, planets on nearby orbits of 
masses ^ M® are expected to produce the largest effect. 
Thus the detection of apsidal precession could actually be 
used to infer the presenc e of Ea rth-mass planets. 

As IJordan fc Bakosl (|2008T ) noted, apsidal precession 
should cause a change in the transit duration and in order 
to estimate the magnitude of this effect, the authors differ- 
entiated Ttso5 with respect to w. Having shown the Ti offers 
substantial improvement over the Ttsos in the previous sec- 
tion, we are here able to provide an improved estimate for 
the secular TDV caused by apsidal precession: 

dT P eo?cosw / 1 . / y/l - b 2 \ 

ir- = — ; ———7— — ^=^= — — / arcsm 

doj 7T (1 - e 2 yi 2 \ y/l^W^a 2 R g 2 - 1 V clrQc sin i ) 

(54) 

We can see that there are two terms counter-acting in 
our derived quantity. The two terms can be understood to 
originate from the planet-star separation changing as a re- 
sult of the precession which has two effects 1) decreasing the 
planet-star separation causes a near-grazing transit's impact 
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parameter to decrease and thus increases T (the first term) 
2) decreasing the planet-star separation causes the tangen- 
tial orbital velocity to increase and thus decreases T (the 
second term). The cos cj term outside of the brackets deter- 
mines the sign of which term causes an increase and which 

to d ecrease . 

iKopall (|l959h showed that t he two effects approximately 
cancel out for b ~ l/y/2. The IKopall l|l959h derivation is 
quite different for the ones produced in this paper. It is done 
by first solving for the mid-eclipse time by a series expansion 
of the differential of the planet-star separation with respect 
to true anomaly, disregarding terms in sin 3 i or higher. The 
duration between the primary and secondary occultation is 
then solved for in another series expansion in first order of 
e. Nevertheless, setting b to this value, the terms inside the 
bracket of equation (54) become: 

= = — 2 arcsin ( = ] (55) 

Under the condition a> JJ,, we find that equ ation (54) 
gives dT/duj = 0, in agreement with lKopall (Il959l ). For very 
close-in orbits, this does not appear to hold. 

We may com pare our estimate of the apsidal precession 
to equation (15) of ljordan fc B akos (2008) , which was found 
by differentiating the expression of TS05 with respect to u). 
The difference between the two expression is typically less 
than 1% across a broad range of parameters. However, if 
b ~ l/v2, the difference between the two diverges and can 
reach 10%- 100%. Given the sensitivity of both equations to 
this critical value of b, we recommend numerical calculations 
over analytic approximations if b is known to be close to 
0.707. 



7.3 Falling Exoplanets 

Planetary bodies experience infall towards the host star 
through tidal dissipation and to a much lesser degree gravi- 
tational radiation. The effects increase as the orbit becomes 
smaller leading to runaway fall-in. Therefore, for very close- 
in exoplanets, the change in semi-major axis may be de- 
tectable. The transit duration will vary as: 

dT P Q 2 C ( 3 . / yi^F \ 

— — ^=^= — arcsm — 

oa 7r ay/i - e 2 \ 2 \a R g c sim/ 



As for apsidal precession, there are two countering com- 
ponents which are the same as before except for a slightly 
different constant in front of the first term. This different 
constant means that the impact parameter at which both 
effects cancel has now changed to b ~ 1/^3 = 0.577. This 
result could not be found in the previous literature and is of 
particular interest given the recent discovery of exoplanets 
on periods of around a day or less, for example WASP-18b 
(Hellier et al. 2009) with period of 0.94 days and b = 0.25. 



7.4 Eccentricity Variation 

Irregular satellites are known to exchange orbital inclina- 
tion and eccentricity through the Kozai mechanism, which 
roughly conserves the value cos Iy/1 — e 2 , where I is the an- 
gle to the ecliptic. Changes in orbital eccentricity are pre- 
dic ted to lead to lo ng-term transit time variations (L-TTV) 
bv iKippinsl (|2008l ). but here we consider the effect on the 
transit duration too. 



7.2 Nodal Precession 

Nodal precession causes changes in the orbital inclination of 
the planetary orbit, which would be a source of s ecular TDV. 
The secular theory o f lMurrav fc DermottJ jl999h predicts the 
rate of inclination change due to a perturbing planet as the 
nodes precess: 



di dui 



(56) 



Where Af2 s k y is the ascending node of the perturb- 
ing planet relative to the ascending node of the transiting 
planet, measured clockwise on the plane of the sky. Thus 
any occurence of apsidal precession from a perturbing planet 
will, in general, be coupled with nodal precession. We may 
derive the rate of secular TDV from inclination change as 
before and find: 



dT 
~di 



- 1 



vr taniVl — e 2 \Jl — b 2 



(57) 



This expression has only one term and therefore we can 
see that decreasing the inclination towards a more grazing 
transit always yields a shorter transit duration, and vise 
versa. 



dT 
~de 



TV (1 



2)5/2 



2e + (1 + e 2 )sinw 



[3e + (2 + e 2 ) sin a;] arcsin 



a,RQ c sini 



(59) 



The two terms here seem to exhibit a more complicated 
inter-dependency which is physically based on the same idea 
of varying the planet-star separation. The balance-point be- 
tween the two effects occurs for: 



e + sin vj 



(3e + (2 - 



(60) 



8 CONCLUSIONS 

We have derived and tested a new approximate expression 
for the transit duration of an extrasolar planet with non- 
zero orbital eccentricity (equation 15). The expression has 
been shown to analytically reduce down to the exact expres- 
sions for a circular orbit, unlike the most previously utilized 
equation. In numerical tests, the new equation is shown to be 
more accurate than the other candidate expressions consid- 
ered in this work, in particular for highly eccentric systems. 



Investigations of transit duration expressions 13 



Quantitvely, the new expression yields a > 200% improve- 
ment in accuracy over the previously most utilized expres- 
sion. 

Manipulation of the new expression provides for a new 
lightcurve fitting parameter set which i s based upon a mod - 
ification of a previously proposed set bv lBakos et al.l l|2007h . 
The new parameter set is shown to demonstrate the low- 
est mutual correlations compared to other most commonly 
used parameter sets and therefore yields the most efficient 
algorithm for fitting lightcurves. 

Additionally, we have shown that the effect of even mild 
eccentricity can cause very large biases in the lightcurve de- 
rived stellar density, which is often used a method for dis- 
criminating between planets and blends in transit surveys. 
Consequently, planetary candidates can be either falsely re- 
jected or accepted for systems with poor contraints on the 
eccentricity. 

Finally, we have used the new equation to derive the 
rates of secular transit duration variation (TDV) as a re- 
sult of apsidal precession, nodal precession (e.g. due to a 
perturbing planet), in- falling extrasolar planets and eccen- 
tricity variation (e.g. Kozai mechanism). These derivatives 
will provide for a more accurate interpretation of secular 
TDV. 
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14 D. M. Kipping 

Table 1. List of important parameters used in this paper. 



Parameter Name Defintion 



5 Sky-projected separation Sky-projected separation of the companion's centre 

and the host star's centre in units of stellar radii 

R„ Radius of the host star Radius of the host star 

Rp Radius of the companion Radius of the companion 

p Ratio-of-radii Ratio of the companion's radius to the stellar radius (Rp/R t ) 

<5 Geometric depth The observed transit depth in the absence of limb darkening and 

blended contamination, equal to p 2 . 

tj First contact Instant when 5 = 1 + p and dS/dt < 

tn Second contact Instant when 5 = 1 — p and dS/dt < 

tc Mid-transit time Instant when dS/dt = i.e. inferior conjunction 

tin Third contact Instant when 5 = 1 — p and dS/dt > 

tiv Fourth contact Instant when 5 = 1 + p and dS/dt > 

tp Total duration Time for companion to move between contact points I and IV 

tp Full duration Time for companion to move between contact points II and III 

T Transit duration Time for companion to move across the stellar disc with entry 

and exit conditions defined as S = 1 

W Transit width Mean value of tp and tp 

tl2 Ingress duration Time for companion to move between contact points 1 and 2 

*34 Egress duration Time for companion to move between contact points 3 and 4 

t Ingress/Egress duration For circular orbits, t\i = <34 = r 

Ti Ti duration A one-term expression for T derived in this work 

T2 T2 duration A two-term expression for T derived in this work 

Ttsos T S05 duration Expression for T derived bv Tinglev fc Sackett (2005) 

T\yio W10 duration Expression for T derived bv Winn (2010) 

^SMO03 SMO03 duration Expression for T derived by Scagcr & Mallcn-Ornclas (2003) 

/ True anomaly True anomaly of the companion during its orbit around the host star 

E Eccentric anomaly Eccentric anomaly of the companion during its orbit around the host star 

M Mean anomaly Mean anomaly of the companion during its orbit around the host star 

H Reduced mass Reduced mass of the companion-star system 

J Angular momentum Angular momentum of the companion 

a Semi-major axis Semi-major axis of the companion's orbit 

ap Semi-major axis Semi-major axis of the companion's orbit in units of stellar radii 

e Eccentricity Orbital eccentricity of the companion's orbit 

f2 Longitude of the ascending node Longitude of the ascending node of the companion's orbit 

to Longitude of pericentre Longitude of pericentre of the companion's orbit 

ui Argument of pericentre Argument of pericentre of the companion's orbit (oj = to — Q) 

q Companion-star separation Companion-star separation in units of stellar radii 

q c Mid companion-star separation Companion-star separation in units of stellar radii at the moment of mid-transit 

P Period Orbital period of the companion 

b Impact parameter Value of 5 when dS/dt = 

i Inclination Orbital inclination of the companion's orbit relative to the line-of-sight of the observer 

D Duration function A parameter defined by Kipping (2008) 

f a True anomaly 'a' The true anomaly of the companion when 5 = 1 and dS/dt < 

f b True anomaly 'b' The true anomaly of the companion when 5 = 1 and d5/dt > 

f c Mid-transit true anomaly The true anomaly of the companion when d5/dt = 

/m Mean transit true anomaly Mean of and f a 

fll Half true anomalistic duration One half of the difference between /& and f a 

E a Eccentric anomaly 'a' The eccentric anomaly of the companion when 5 = 1 and dS/dt < 

Eh Eccentric anomaly 'b' The eccentric anomaly of the companion when 5 = 1 and dS/dt > 

Ejj Mean transit eccentric anomaly Mean of E b and E a 

En Half eccentric anomalistic duration One half of the difference between E b and E a 

p t Stellar density Average density of the host star 

C/R* Zeta over R„ The reciprocal of one half of the Ttsos duration. 

T / R* Upsilon over R * The reciprocal of one half of the Ti duration 

T> Deviation Deviation of a candidate duration expression from the exact solution 

I Improvement Improvement of a candidate duration expression over the Ttsos expression 



